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Abstract 



^Sj . We present a separable expansion approximation method for Coulomb- 

like potentials which is based on Schwinger variational principle and uses 

■^-j- \ Coulomb-Sturmian functions as basis states. The new scheme provides faster 

convergence with respect to our formerly used non-variational approach. 



Both variational approaches and separable expansion schemes are extensively used in 
solving few-body problems. Some time ago Adhikari and Tomio presented an unified treat- 
ment of separable expansion schemes based on Schwinger variational principles |T||. They pro- 
posed various approximation schemes for the transition operator t satisfying the Lippmann- 
Schwinger equation 

t = v + vg°t, (1) 



^ . where v is the potential and g is the free Green's operator. It was found that using 
these schemes with appropriate choice of expansion functions a rapid convergence could be 
obtained. However, in this paper, not a single word was devoted to Coulomb-like potentials. 

At about the same time in a series of papers a separable expansion scheme for 
Coulomb-like potentials was proposed by one of us. The Coulomb interaction was kept 
in the Green's operator and only the short range part of the potential was subject to the 
separable expansion. This approach uses Coulomb-Sturmian (CS) functions as basis allow- 
ing an exact analytical calculation of the matrix elements of the Coulomb Green's operator 
(see ||, and recently ||). Thereby only the short-range part of the interaction is approx- 
imated and the correct Coulomb asymptotic properties of all quantities are guaranteed. 
The method has also been applied in three-body Faddeev calculations for bound-state and 
scattering problems with Coulomb interactions M. 

In this paper we generalize one of the separable expansion schemes proposed by Adhikari 
and Tomio for two-body Coulomb-scattering problems. 

The expansion schemes in Ref. |IJ are based on a finite rank N approximation of the 
product of operators 
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N 

AB- 1 C^Y, A h)D lJ (QC, (2) 

where 

{D- 1 ) ji ={( J \B\r H ). (3) 

Here the sets \r]i) and \Q) are assumed to be complete sets of states in the Hilbert space. If 
N tends to infinity Eq. (Q) becomes an identity. 

From among the separable approximation schemes of Ref. |J let us start with the form 
of t following from Eq. ([J) 



t=l{v~ 1 -go)' 1 !. (4) 



Applying Eq. (0) we get 



N 



with 

(P~ X )ji = (Vjlv^h) ~ (vMVi)- (6) 
The first term in Eq. (^|) can again be approximated as before 

(r)j\v~ l \r)i) = (r] j \lv- 1 l\r] i ) 

N 

(7) 



where 



(c- 1 )^ = (d'\v\c f ). (8) 



This approximation scheme can easily be generalized for Coulomb-like potentials. A 
Coulomb-like potential in some partial wave I can be written in the form 

vi = v c + vf, (9) 

where v is the pure Coulomb potential and vf is short ranged. The Lippmann-Schwinger 
equation for the Coulomb-modified transition operator t sC reads 

tf = v! + v?gf(E)tf, (10) 

where gf(E) = (E — h® — v )' 1 is the Coulomb Green's operator and h® is the free Hamil- 
tonian. The solution of this equation can be given in a form analogous to Eq. (H) 



tf = l 



(vtr'-gfiE) i, (ii) 
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and the whole procedure of Eqs. ©-© can be repeated. Only the free Green's operator 
has to be replaced by the Coulomb one. 

As basis states we choose CS functions because they allow an exact analytical calculation 
of matrix elements of gf(E). In coordinate representation they have the form 



(r|n; b) 



(n + 2l + l)\ 



1/2 



(2br) l+1 exp (-br)L% +1 (2br 



(12) 



where n = 0,1,2,..., L^+i 

are the Laguerre polynomials and b is a scaling parameter. They 
have an analogous simple form also in momentum representation 



(p\n; b) 



2 l +VH\{n + I + 1] 



n(n + 2l + l)\ 



fe(2M m rl+1 ( P 2 ~b 2 ' 



(13) 



where G denotes the Gegenbauer polynomials. Introducing the notation (r|n; b) = (r\n; b)/r 
we can express the orthogonality and completeness of the CS functions as 



(n; b\n'; b) = 5 n 



and 



N 



1 = Jim ^ \n; b)(n; b\ = Jim ^ |n; b)(n; b\ 



n=0 



n=0 



respectively. 

As \r]i) and \Q) we take CS bases with different b parameters: 



tf ~ Yl |w;&i)A» n /(n';&i| 

n,n'=0 



(14) 



(15) 



(16) 



N 



(D 1 ) n 'n= (n';b 1 \m;b2}C mm i(m';b 2 \n;b 1 ) 

m t m'=0 



■{n']b 1 \g l |n;6i) 



(17) 



(C- 1 ) m >m = (m';b 2 \vf\m;b 2 ). (18) 

While the matrix elements of the potential and the overlap of CS functions should be cal- 
culated numerically either in configuration or momentum space the matrix elements of the 
Coulomb Green's operator can be calculated analytically |2|||. 

This scheme is essentially equivalent with a twofold separable expansion for potential vf 
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N 



~ Yl \n;b 1 )((n;b 1 \m';b 2 )) 1 



n,n' ,m,m'=0 

(m; b 2 \vf\m; b 2 ) 

{(m^b^n'-h))- 1 ^^. (19) 

It is easy to verify that Eqs. QlBQ- fllSD are the solutions of the Lippmann-Schwinger equation 
Eq. ([TOD , with this approximate potential. 

This scheme should be compared with the CS potential separable expansion scheme of 
Refs. Pfl which is not supported by variational principle. It is based on the approximation 
of the unit operator 



1= lim £l^nH (20) 

n=0 

with o factors possessing the properties lim^oo = and limjv^oo^ = 1- Now, the 
approximation takes the form 



N 



V 



lvflw £ \n)a»{n\v sl \n')a»{n>\. (21) 



n,n'=0 



For the form 



N l-exp{-[a(n-iV-l)/(iV + l)] 2 } 



l-exp(-a 2 ) 1 J 

was used with a ~ 5. This value of the arbitrary parameter a yielded the fastest convergence. 

To show the relative power of these separable expansion schemes we have calculated the 
Coulomb-modified nuclear phase shifts 5f c at / = and at various energies for a p — p 
scattering. The short range potential was taken in Malfliet-Tjon form 

Vi = v exp(-f3 r)/r + v 1 exp(— /?i?-)/r, (23) 

with v = -626.885MeV, (3 = 1.55fm, v x = 1438.720MeV, fa = 3.11fm. In the method 
of Eq. (|2~ID CS basis parameter b = 3fm _1 was used, while in the expansion Eq. ( |19|) bi = 
3.8fm _1 and b 2 = 2.5fm _1 were taken. It can be seen in Table I that in both approximation 
schemes it is possible to choose the b parameters so that the expansions give almost equally 
fast convergence over the whole spectrum and provide extremely accurate results. To reach 



6-digits accuracy the method of Eq. (|1^) needs 10-13 basis states, while the method of Eq. 
(pl|) needs 20-23 states. We have observed similar results over a wide range of b parameters. 
It should be noted, however, that the method of Eq. ([19|) is more complicated numerically, 
so the numerical effectivity of both methods are more or less the same. 

In this Brief Report we have combined two separable expansion methods. One based 
on Schwinger variational principle has been proposed by Adhikari and Tomio in Ref. [|TJ. 
The other approach proposed by Papp in Refs. is not variational and it was designed for 
Coulomb-like potentials. This new scheme is a variational separable expansion method and 
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is applicable for Coulomb- like potentials. It converges considerably faster in terms of basis 
states then the non- variational method of Eq. (PH) . This property could be useful in three- 
body calculations where the rank of the expansion is of crucial importance. The method 
of Eq. (pip has been generalized for solving Faddeev-type integral equations of three-body 
Coulombic systems. Whether or not the method of Eq. ([0|) can be extended in this direction 
is still an open question. 

This work has been supported by the OTKA Gratns No. T026233 and T029003 and the 
NSF Grant No.Phy-0088936. 
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TABLES 



N E = 


0.1 MeV 


E = 


1 MeV 


£ = 


10 MeV 


E = 100 MeV 


Eq. © 


Eq. (IS) 


Eq. (|21|) 


Eq. (|19|) 


Eq. © 


Eq. © 


Eq. (|l]) Eq. (|l|) 



2 


-0 


100682 


-0 


151364 


-0 


624200 


-0 


797903 


1 


546662 


1 


453638 





163942 





253576 


3 


-0 


115923 


-0 


122391 


-0 


685037 


-0 


710521 


1 


506720 


1 


479076 





365679 





351795 


4 


-0 


120857 


-0 


120468 


-0 


706248 


-0 


705228 


1 


480305 


1 


477592 





430684 





403522 


5 


-0 


118844 


-0 


119136 


-0 


701335 


-0 


701377 


1 


473793 


1 


477278 





415215 





402637 


6 


-0 


118655 


-0 


119182 


-0 


699997 


-0 


701573 


1 


481240 


1 


479584 





411720 





404552 


7 


-0 


119150 


-0 


119165 


-0 


701393 


-0 


701565 


1 


483505 


1 


480309 





408933 





407128 


8 


-0 


119166 


-0 


119165 


-0 


701534 


-0 


701549 


1 


481915 


1 


480743 





408615 





407224 


9 


-0 


119147 


-0 


119164 


-0 


701390 


-0 


701523 


1 


481590 


1 


480921 





409217 





407345 


10 


-0 


119155 


-0 


119164 


-0 


701439 


-0 


701519 


1 


481230 


1 


480945 





408458 





407467 


11 


-0 


119157 


-0 


119162 


-0 


701462 


-0 


701508 


1 


481028 


1 


480958 





407939 





407487 


12 


-0 


119158 


-0 


119162 


-0 


701472 


-0 


701507 


1 


480981 


1 


480957 





407830 





407488 


13 


-0 


119160 


-0 


119162 


-0 


701490 


-0 


701505 


1 


480959 


1 


480957 





407694 





407497 


14 


-0 


119160 


-0 


119162 


-0 


701492 


-0 


701504 


1 


480958 


1 


480957 





407585 





407499 


15 


-0 


119161 


-0 


119162 


-0 


701498 


-0 


701504 


1 


480959 


1 


480957 





407557 





407499 


16 


-0 


119161 


-0 


119162 


-0 


701499 


-0 


701504 


1 


480959 


1 


480957 





407535 





407499 


17 


-0 


119161 


-0 


119162 


-0 


701501 


-0 


701504 


1 


480959 


1 


480957 





407511 





407499 


18 


-0 


119161 


-0 


119162 


-0 


701502 


-0 


701504 


1 


480959 


1 


480957 





407506 





407499 


19 


-0 


119161 


-0 


119162 


-0 


701502 


-0 


701504 


1 


480958 


1 


480957 





407506 





407499 


20 


-0 


119162 


-0 


119162 


-0 


701503 


-0 


701504 


1 


480957 


1 


480957 





407500 





407499 


21 


-0 


119162 


-0 


119162 


-0 


701503 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


22 


-0 


119162 


-0 


119162 


-0 


701503 


-0 


701504 


1 


480957 


1 


480957 





407501 





407499 


23 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


24 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


25 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


26 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


27 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


28 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 


29 


-0 


119162 


-0 


119162 


-0 


701504 


-0 


701504 


1 


480957 


1 


480957 





407499 





407499 



TABLE I. Convergence of the Coulomb-modified nuclear phase shift 5q C (E) (in radians) in 
the (^) potential at different energies with respect to the ./V number of basis states used in the 
separable expansion schemes (|2l| ) and (|i"9|). 
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